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Abstract 



We describe flux tubes and their interactions in a low energy sigma model 
induced by SU(Np) — » SO(Np) flavor symmetry breaking in SO(N c ) QCD. 
Gauge confinement manifests itself in the low energy theory through flux tube 
interactions with unscreened sources. The flux tubes mediating confinement 
act as Alice strings in their cores, a phenomenon which may occur for 7T2-line 
defects in physically realized systems. 
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Nonlinear sigma models successfully describe many low energy QCD phenomena. How- 
ever, they have not captured the hallmark feature of QCD: confinement, where potentials 
between qq pairs grow linearly with separation. Such linear potentials arise from flux tubes 
between unscreened q and q sources for pure Yang-Mills QCD. Conventional sigma mod- 
els describe the low-energy dynamics of QCD due to the global flavor symmetry breaking 
SU(Nf)l x SU(Nf)r — > SU(NF)diag- Global flux tubes do not occur in these models, as 
the vacuum manifold G/H forms a Lie group, with trivial ir 2 . This absence of flux tubes is 
consistent with the absence of unscreened sources in QCD with fundamental quarks. 

Witten noted that an SO(N c ) gauge theory of QCD induces a different sigma model, 
whose topology can support flux tubes. The theory has Np left-handed quarks which trans- 
form as (real) fundamentals of SO(N c ). The color singlet condensate (qf^ q^j) (where a and 
i are color and flavor indices respectively) breaks the SU(Np) flavor symmetry to its SO(Np) 
subgroup. There is a baryon number SO(N c ) 2 anomaly which leaves an anomaly free con- 
served Z 2 baryon number. The Goldstone modes are described by a SU(Np)/ SO(Np) non- 
linear sigma model, with skyrmions (ir 3 (G/H) = Z 2 , for Np > 4) and flux tubes {ir 2 {G/H) 
= Z 2 , for Np > 3). This defect classification makes sense: (a) baryons are identified with an- 
tibaryons since real quarks are identified with antiquarks; (b) there are Z 2 flux tubes which 
confine spinor sources, since only an even number of spinors can be screened by fundamental 
quarks. 

In this letter we construct the flux tubes in this theory and show that their interactions 
with skyrmions and spinor sources obey heuristic expectations. We show that confinement 
in an SO(N c ) gauge theory can manifest itself in the low energy sigma model through flux 
tube interactions with unscreened sources. We proceed as follows: we derive the unique 
flux tube form with minimal energy; examine its classical stability and dynamics; and then 
study its quantum stability, spectrum, and interactions. We discuss how flux tubes act as 
Alice strings in their cores, form loops that support skyrmion number, and mediate the 
confinement of spinor sources. We skip many technical details, which appear in |fj|. 

Finding Nontrivial Flux Tubes 
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The sigma model for SO(N c ) QCD describes the dynamics of a sigma field E, which 
encodes the orientation of the fermion condensate (ggj qfj) with respect to its standard form 
oc 5ij, with £ = 1. Under an SU(N F ) transformation by g, E transforms as E — > <?£g T , 
so that a generic £ has the form gg T for some group element g. E = 11 is invariant 
under a g transformation if g E SO(N F ), so E is an element of the coset space G/H = 
SU(Np)/ SO(Np), which is not a Lie group. 

G invariance and cylindrical symmetry fixes the form of the minimal energy flux tube to 
be E = g(9)T> (r)g T (9) — an r-dependent vev, with angle-dependent group rotation, where 
(r,0,z) are cylindrical coordinates. Choosing E c = 11 at spatial infinity forces g{6) to lie in 
H = SO(Np), and to commute with E G (r = 0), for a nonsingular solution with finite energy. 

E can be further specified using basis generators of the su(N F ) Lie Algebra, with Cartan 
norm tr T a T b = \5 a b- The rank N F — 1 Cartan subalgebra has as a basis diagonal matrices 
T d , 

T d = (2d(d + 1) Y 1 ' 2 diag(l, . . . , 1, -d, 0, . . . , 0), (1) 

with ones in the first d entries and 1 < d < N F — 1. We can overspecify them, sacrificing the 
Cartan norm, by Pauli matrices {\r z (jk) } in all (jk) subplanes of the N F dimensional vector 
space. Off-diagonal generators {_\r x (jk)i \ T y{jk)} complete the basis. The antisymmetric 
matrices {T h } = {^T y (jk)} generate the unbroken H = SO(N F ) symmetry, while the sym- 
metric matrices {T b } = {\t x ^) , ^T z Qk) } are broken generators. E G (r), a unitary symmetric 
matrix, can be written as exp {iFb(r)Tb}. Thus E assumes the most general nonsingular 
form E = h{6) exp {iF b (r)T b )}h^ 1 (9), for h{6) G H and F b (r) ranging from zero at infinity 
to 2nn 5bb> (for some direction b') at the origin. 

We construct a non-trivial flux tube E from the exact sequence 

7T 2 (SU(N F )) = -> 7T 2 (SU(N F )/SO(N F )) -»• 7Ti (SO(N F )) = Z 2 -> m (SU(N F )) = 0. 

That is, gg T gives a nontrivial E only if g corresponds to some mapping from the plane to 
SU(N F ), with boundary values in the SO(N F ) subgroup. Furthermore, when parametrized 
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as a family of loops, g must start at the identity and end on a nontrivial loop in SO(Np)- 
Taking (a G [0, 2tt], /3 G [0, 7r]) as our coordinates on the plane, these criteria become 

!1 when a = 0, a = 2n, or /5 = 
(2) 
h 2 (a) when (3 = n, 

where /i 2 (a) is a nontrivial loop in SO(Np). 

Through technical arguments, two of us show that such trivializations of SO(Np) loops 
in SU(Np) have minimal energy only when they induce £ of a very limited form 0. To 
minimize energy, h 2 (a) is geodesic: h(a) = exp {ia T hn }, where T hn = {\r y } in some plane 
□ — taken here as (12) for concreteness. This geodesic loop determines the associated flux 
tube E: 

£(r,0) = h{9) b{r) hr l (B), with 

h{6) = exp {i9 Thn}, b(r) = exp {iF d {r)T d }, (3) 

and Fi(r = 0) = 2tt; F d (r -> oo) = 0. 

Here T/j D generates the loop /z 2 (a), and are the Cartan generators from eq. ([!]). The 
boundary conditions on F d stem from eq. (0), requiring g(a,/3) to trivialize /i 2 (a). 

Note that when F d >\ vanishes, £ lies entirely within a planar SU(2) subgroup of SU(Np), 
as Ti = r z n /2. This produces the simple form 

S n (r, 0) = 11 + ( cos(F/2) - 1 ) 11 D + i sin (F/2) ( cos 6 r zU - sin 6 r xU ) . (4) 

where F = F 1; 11 □ gives the identity in the plane □ and vanishes outside it, and 1 is the 
usual SU (Np) identity. In this letter, we consider only such planar flux tubes, having shown 
in that they minimize energy for the action discussed here. 

Flux Tube Stability and Dynamics 
As a minimal model for S with stable skyrmions, consider the Skyrme lagrangian || [] 

*We use QCD sigma model parameters e,F n , m w , etc., although these are not the usual pions. 
There is also a Wess-Zumino term, which vanishes for constant-z flux tubes. 
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Co = ^ tr 9 M E^Et + J-j tr [E^E, E^E] 2 . (5) 

This model does not have stable flux tubes, for under the rescaling E(r, 9, z) — ► S(Ar, 0, z), 
the tension of a finite flux tube can decrease. Flux tubes diffuse to infinite size (A = 0) to 
lower their energy. 

There are several modifications to the minimal Cq, or to the minimal form for E, that 
lead to stable flux tubes 0. One natural choice is to give non-zero mass to the quarks in 
the original QCD theory. This induces a mass term in the effective lagrangian, 

F 2 m 2 

£ m = ^^tr (E + Et-2-ll) , (6) 

which gives mass to the Goldstone bosons, and also stabilizes the flux tube solutions. 
The Skyrme action eq. (|5]) gives gradient energy density to the flux tube: 

Po = §tr {(i^) 2 + l (T 2 -[^T d ,T] 2 )} , (7) 

where 

f = b- l {r)[T m , b(r)}, (8) 

and r has been rescaled to dimensionless units eF^rphys- For the planar solution eq. (f|), 
this gives the energy density 

Po = f {^ /2 + l(l-cosF)(l+F' 2 )} . (9) 

To this term we add potential energy from the pion mass term eq. (Q), 

Pm = ^ (1 - cos(F/2)) , (10) 

where A = 2m. n /e F n . Together, these energy contributions determine a nonlinear equation 
of motion for F, which we solve numerically. Solutions appear in Figure 1 for different 
values of A, including the physical Ao = 0.236 (e = 2tt, m n = 138 Mev and F n = 186 Mev). 
Increasing A raises the flux tube's energy density while shrinking its core size. Inside the 
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core, F falls linearly from F(0) = 2tc; outside, it scales as the hyperbolic Bessel function 
F ~ x~ x l 2 exp{— x}, with x = Xr/2. 2A _1 sets the flux tube's dimensionless core size, 
which gives physical core size m~ x . The flux tube tension is proportional to F% — which is 
proportional to the number of colors N c in the large N c limit. This supports its identification 
with the confining force between spinor sources, since the Casimir for SO(N c ) spinors is of 
order N c . Numerically we find tension 4.6F% when A = Ao, scaling linearly with A. 

Quantum Stability and Spectrum 
The quantum numbers and low-lying internal excitations of the flux tube are given by 
quantizing the zero-mode collective coordinates @, 

E(t,r,0)=A(t)E(r,9) A~\t) , (11) 

where A(t) rotates about n^T^ with dimensionless frequency u. Q These modes have rota- 
tional energy confined to the string core, which can be calculated from the Skyrme action 
eq. (§: 

A- = ^tr (T 2 -[FiT d ,f] 2 -l[T,f] 2 ) . (12) 

with T from eq. (§) and 

f = b-\r) [ h-\6) n h T h h{6) , b(r) ] . (13) 

Calculating eq. (]T2|) is tedious, and described in 0. We show there that the planar vacuum 
survives quantum fluctuations due to zero modes. It has the classical two dimensional 
Lagrange density 

L(z,t) = -p - ^(tr A^AT hD f - ^^(tr ^ AT ha ,)\ (14) 

where z measures dimensionless length along the flux tube and the planes □' intersect □ 
in single lines. Our numerical solutions, modified by rotational back reaction, give for the 



^We neglect other excitations, such as bending modes. 
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integrals p ~ 5F n /e, A n ps 40/e 3 F 7r , and Ap/ ~ 16/e 3 F.,,. The moments of inertia A n and 
A a / vary little with back reaction for quantum excited states; however, the tension p grows 
linearly with u, due to compression of the rotating flux tube. 



To quantize this Lagrangian, we must rewrite the time derivatives in eq. (|T^) in terms 
of canonical momenta. Under an SO(Np) transformation h, £ — > hHh T , so that A — > hA 
from eq. (0). The generators 1h of H rotations thus generate left-transformations of A. 
Similarly, we define body-centered generators T' h generating right transformations, A — > A/i. 
Xft and 1' h are related by an orthogonal transformation, 

l h = R hh ,l' h ,, R hh , = 2 tr (A*T H AT h ) . (15) 

Noether's theorem then gives for X' h , 

l' h = iA h tr {A^ATh) (no sum). (16) 

The flux tube solution E(r, 9) is invariant under the action of T' D + J z , so that T' D is — J z 
(much like the Skyrme model, where / in the body centered frame is — J). It is also invariant 
under SO(N F — 2) C SO(N F ), acting on the subspace Dj_ orthogonal to □. This restricts 
S to quantum states X' containing a singlet under SO(Np — 2). 

Given these constraints, we write the Hamiltonian density obtained from eq. (|14|) in 
terms of the physically observable Noether charges, and J, 

mt) = p + ^j 2 z + ^s i2 - J t ) . ( l7 ) 

which is quantized subject to the constraint 1' D = —J z , with bosonic defects for Np > 3. 
For three flavors, this gives a spectrum of states (J, J z ) where / and \ J Z \ are integers, —I < 
J z < /, and J 2 = 1(1 + 1). 

Flux Tube Interactions 

Our sigma model has both flux tube and skyrmion solutions. We show that flux tube 
solutions are Alice strings. Moreover, loops of flux tube can carry skyrmion number, just 
as loops of gauged Alice string support monopole charge. We construct the skyrmions 



and discuss their interactions with flux tubes, showing that only the topologically trivial 
combination of two flux tubes can end on a skyrmion. This suggests that, while baryons 
(i.e. skyrmions) are not confined, the spinor sources which combine to form them are, with 
confinement mediated by Z2 flux tubes joining them. 

First we consider twisted flux tubes. Our flux tubes have the Alice property: some 
unbroken symmetries preserve a local vev but cannot be extended globally, since transport 
around the string makes them multivalued Q. This property seems ambiguous, as the flux 
tube eq. (||) allows many choices of underlying g(r,9). Double- valued choices — such as 
g(r,6) = h(9) b(F(r)/2) — give unbroken generators g Thw g^ 1 double-valued in 9; while 
single-valued choices — such as g(r,6) = h{9) b(F(r)/2) h{9) — leave all generators single- 
valued. This ambiguity persists even when we recast the flavor-dependent quark mass S 
as an interaction between flavor gauge fields g~ x {r : 9) <9 M g(r,9) and shifted quarks Ql = 
g(r,9) qL. || This identifies the flux tube with a gauged string, with Wilson loop U(2ti) = 
Pexp [§ A - dl However, the ambiguity persists: g(r,9) double- valued gives Wilson loop 
U(2ti) = 1 — 2 11 □, making both quarks Ql and generators T hn i double- valued. Single- 
valued g(r, 9) produces instead single- valued quarks with an r-dependent Wilson loop, giving 
U(2ti) = 1 and single- valued generators at r = 00, but nontrivial U{2i{) and multi-valued 
generators Thw at finite r. 

This ambiguity can be resolved physically, by considering adiabatic transport of quarks 
around the flux tube. Under such transport, quarks remain in their mass eigenstates. 
These are fixed by two terms: a flavor-independent bare Majorana mass M (breaking 
SU(Nf) — > SO(Nf) explicitly) and the flavor- dependent Majorana mass /zE, where /i is 
the S , [/(A^/)-breaking vev. These give mass eigenstates double-valued in 9. However, the 
states have mass splitting Am 2 = —2iM(fi — fi*) sin(F/2), with M real. Thus the quarks 
are degenerate, and unaffected by transport around the string, unless the vev li misaligns 
in phase with the bare mass M. In that case, quarks at finite radius have double-valued 
wave functions and Aharonov-Bohm scatter off the flux tube, flipping the sign of their T^w 



S 



charges. Thus the flux tube acts as an Alice string in its core. This analysis readily general- 
izes to other 7r 2 -strings. This suggests the promise of observing similar effects in nonabelian 
physical systems with global 7r 2 -strings, such as liquid He-3a. J7J 

Since twisted Alice loops can support monopole charge, we calculate the skyrmion num- 
ber of a twisted flux loop, £(z, r, 9) = A(z) E(r, 0) A^ x {z) with A(z) = exp [iz I n^T/j) . 
27r-periodicity implies that I is integral for planar = n n , and even otherwise. Thus any 
twisted flux loop can deform to the planar flux loop £(z, r, 6) = h(9 + lz) b(F(r)) h~ 1 (9 + lz), 
taking values within the subspace SU{2) / SO{2) ~ S 2 . Thus we can identify its 7r 3 index 
with the map's Hopf number, i.e. the linking number between any two fibers of constant 
E in physical space. From ||, this linking number is precisely I — the number of times 
a nontrivial fiber S twists around the loop's core, which has S = — 11 □. Thus flux loops 
with an I = 1 planar twist form fundamental skyrmions; flux loops with nonplanar twist are 
trivial for Nf > 3. 

A nicer parametrization of the skyrmion stems from the exact sequence 

7r 3 (SO(N f )) - 7r 3 (SU(N f )) - 7r 3 (SU(N f )/SO(N f )) - 7r 2 (SO(Nf)) — . 

This identifies skyrmions on the vacuum manifold with images of skyrmions in SU(Nf), 
g(r,h) = exp(iF s (r) hi Tin)- (Here r and h are the radius and unit direction vector in 
3-space, and F s (r) approaches 2tt at r = and zero at r = oo.) This gives an axisymmetric 
skyrmion 

S s = 11 +(cosF s -l) (l-n 2 z ) lu 

+i (sin F s (n y r xU + n x t zU ) + 2 sin 2 (F s /2) n z {-n x t xU + n y r xa )) , (18) 

after global spatial rotation. Equation (|J) for the flux tube lets us identify the angular 
winding of a flux tube with that of a skyrmion in the xy-plane (n z = 0). Thus, if their 

^We treat quarks with trajectories well-separated from the origin. This holds asymptotically inside 
the flux tube, where m~ l < r < m^ 1 , for quark mass eigenvalues . 
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radial forms coincided, we could deform the skyrmion's lower hemisphere into a flux tube. 
However, since both F and F s vary from to 2ir, the skyrmion cannot end in a single 
flux tube. Instead it joins only to flux tubes where F(r) ranges from to Att — that is, 
configurations with two flux tubes, which are trivial. Thus skyrmions are not confined. 

However, objects which combine to form skyrmions can interact with the flux tubes. 
Such "half- skyrmions" could arise as external spinor sources in the underlying theory. They 
are confined, as fundamentals cannot screen them. As mappings on G/H, they appear 



precisely as half-skyrmions — objects of form (|Tq) with F s (r) ranging from to tt. Such 
objects are not defects in the conventional sense, since they have linearly divergent energy 
- like an unscreened point source. They can join to their opposite winding counterparts 
via single flux tubes, restricting linearly divergent energy to their separation distance. Thus 
confinement of sources in an SO(N c ) gauge theory persists in the low energy Skyrme model. 
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FIGURES 

FIG. 1. a.) Flux tube solutions F(r) for A values 0.2,0.236,1.0, and 2.0. The dotted line 
corresponds to the physical value Ao = 0.236. Core size shrinks with increasing A. b.) The above 
flux tubes' energy density p (po + p m , from text). Tension grows with A. 
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